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Abstract 

A class of models of biological population and communities with a singular equi- 
librium at the origin is analyzed; it is shown that these models can possess a dy- 
namical regime of deterministic extinction, which is crucially important from the 
biological standpoint. This regime corresponds to the presence of a family of ho- 
moclinics to the origin, so-called elliptic sector. The complete analysis of possible 
topological structures in a neighborhood of the origin, as well as asymptotics to 
orbits tending to this point, is given. An algorithmic approach to analyze sys- 
tem behavior with parameter changes is presented. The developed methods and 
algorithm are applied to existing mathematical models of biological systems. In 
particular, we analyze a model of anticancer treatment with oncolytic viruses, a 
parasite-host interaction model, and a model of Chagas' disease. 

Keywords: Non-analytic equilibrium; ratio-dependent response; pathogen trans- 
mission; elliptic sector; population extinction 



1 Introduction 

Mathematical models of predator-prey or parasite-host interaction and epidemiological 
models formulated as a system of ordinary differential equations (ODEs) share the same 
framework of modeling. The principles of these models, balance relations and decom- 
position of the rates of change into birth and death processes, which were applied since 
publications of Lotka-Volterra equations [H |2] and ODE SIR-model of Kermack and 
McKendrick [3], have remained valid until today. Modifications were limited to replacing 
growth, death, and transmission rates by more complex functions, other types of func- 
tional responses, density- dependent mortality rates, or modes of pathogen transmission 
other than mass action kinetics. 
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In a number of cases different population models possess a distinct similar feature that 
complicates the analysis of the models; namely, a number of models are not defined at 
a singular point. Without loss of generality we can assume that this point is the origin 
0(0, 0). The models we consider in the present paper are formulated in such a way that 
this singularity is removable, for example, we can get rid of it by a time change. After 
the time change a topologically equivalent dynamical system in Rl is obtained, where 
M.2 = {{^jV)'- X ^ 0, y ^ 0}, which is a natural state space for biologically motivated 
models. The origin 0(0, 0) becomes a well defined equilibrium for this dynamical system. 
Sometimes, however, the usual approach by linearization fails to infer the structure of 
a neighborhood of 0(0,0) because the origin is a non-hyperbolic equilibrium for any 
parameter values (i.e., this point has both eigenvalues equal zero). 

The main goal of the present paper is to describe completely the possible structures of 
a small neighborhood Q of non-hyperbolic equilibrium point 0(0, 0) in case of a particular 
class of differential equations that includes many biological models. We also present an 
algorithm to analyze such equations in Q and apply this algorithm to a number of math- 
ematical models. In particular, we show that the models under consideration can possess 
a dynamical regime of deterministic extinction, which is crucially important from the bi- 
ological standpoint. This regime corresponds to the presence of a family of homoclinics 
to the origin, so-called elliptic sector. We note that the fine structure of the phase plane 
containing the elliptic sector was sometimes overlooked in the analysis of these models. 

More specifically, we analyze non-hyperbolic equilibrium 0(0,0) of the following sys- 
tem of ODEs: 

dx 

— = P2{x,y) + P*{x,y), 
dy 

— = Q2{x,y) + Q*{x,y), 

where P2{x,y), Q2{x,y) are homogeneous polynomials of the second order: 

P2{x, y) = p2,ix^ + Pi,2xy + Po,32/^ 
Q2{.x, y) = q^^QX^ + q2,ixy + gl,2?/^ 

and P*{x,y) = 0{\{x,y)\^), Q*{x,y) = 0{\{x,y)\^), i.e., their Taylor series start with 
terms of order three or higher. 

It is important to stress that many of the global properties of the models with the dis- 
cussed peculiarity are determined by the properties of equilibrium 0(0, 0), and, therefore, 
it is essential to know the exact structure of a neighborhood of this point and dependence 
of this structure on the model parameters. 

We organize the paper as follows: Section 2 is devoted to examples of biologically 
motivated models, which, after a suitable time change, fall into the class of differential 
equations ([T]); in Section 3 we introduce the basic definitions and notations and state the 
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main theorems that describe possible topological structures of the origin of system ([T]); 
here we also present an algorithm to study the structure of fl of 0(0,0); in Section 4, 
using the algorithm from Section 3, we analyze some mathematical models, mainly those 
introduced in Section 2; Section 5 contains discussion and conclusions; finally, proofs of 
some mathematical statements are given in Appendix. 

2 Motivation and background 

Ratio-dependent models. Deterministic predator-prey models can be written in the 
following 'canonical' general form: 



where N and P are the densities (or biomasses) of prey and of predators, respectively. 
The production of prey in the absence of predators is described by the function F{N), 
whereas g{N, P) is the functional response (number of prey eaten per predator per unit 
time 1^). The constant e is the trophic efficiency, and predators are assumed to die with 
a constant death rate q. 

Many questions in predator-prey theory revolve around the expression that is used for 
the functional response g{N,P). Much early work was only concerned with the way in 
which this function varies with prey density (e.g., the so-called HoUing types I, II, III), 
ignoring the effect of predator density. For example, in the Lotka-Volterra models, where 
g{N,P) = aN, or in models with the Holling type II response [5j, where g{N,P) = 
aN / {1 + ahN), the functional response is of the form g = g{N) (termed 'prey-dependent' 
by Arditi and Ginzburg [6]), see also [7j. 

It was also recognized that the predator density could have a direct effect on the func- 
tional response. A number of such 'predator-dependent' models have been proposed, the 
most widely known are those of Hassel and Varley ^ , DeAngelis et al. ^ , or Beddington 



Arditi and Ginzburg suggested that the essential properties of predator dependence 
could be rendered by a simple form which was called 'ratio-dependence'. The functional 
response is assumed to depend on the single variable N/ P rather than on the two separate 
variables N and P. Under this supposition system ([2D becomes 



dN 

Itt 
dP 

'dt 



F{N)-g{N,P)P, 



eg{N,P)P-qP, 



(2) 



dN 

It 
dP 

~dt 



F{N) - g{N/P)P, 



eg{N/P)P - qP. 



(3) 
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Discussion of the biological implications and relevance of the ratio-dependent models, 
together with their principal predictions, can be found in, e.g., [TT| \n\. 
One particular model that received considerable attention is of the form 



dN 
"dt 

dP 

HI 




) 



P + ahN ' 



aNP 



(4) 



where Michaelis-Menten or Holling type II function 



g{z) = az/{l + ahz), z = N/P 



was used as the functional response. 

The model Q was studied in, e.g., [H [II [ISl [ISl [HI [IE]. It was shown, by this 
particular example, that ratio-dependent models can display original dynamical properties 
that cannot be observed in two-dimensional prey-dependent models. For instance, the 
origin can be an equilibrium point simultaneously attractive and repelling, thus shedding 
light on ecological extinction. Coexistence of several dynamical regimes with the same 
set of parameters was also observed. 

From the mathematical point of view, ratio-dependent predator-prey models raise 
delicate questions because the functional response is undefined at the origin = 0, P = 0. 
As a consequence, the origin is a so-called non-analytical complicated equilibrium point 



We are only interested in the dynamics of system ^ in = {(A^, P) : A^ ^ 0, P ^ 0}. 
Time scale change dt ^ {P + ahN)dt for system (jl]) results in the system 



The system ([5]) is well defined at 0(0, 0), however, the Jacobian matrix at 0(0, 0) is a zero 
matrix, i.e., 0(0,0) is a non-hyperbolic equilibrium of ([5]). Obviously, model ([5]) belongs 
to the class of ODEs 

Frequency-dependent pathogen transmission in a population of variable size. 

Transmission is the key process in a host-pathogen interaction. In most early models 
transmission was assumed to occur through the law of mass-action. For example, in the 
classical model of two differential equations, formulated as a particular case of the general 
model of Kermack and McKendrick in their pioneer work [3J, the transmission function 
takes the form PSI. Here S is the density of susceptible individuals, / is the density of 



dN 

dP_ 
'dt 



rN{l - N/K){P + ahN) - aNP, 



eaNP-qP{P + ahN). 



(5) 
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infective individuals, and (3 is the transmission coefficient. Note that, for mass action, the 
contacts per unit time per individual rise linear with the population density N = S + I 
(careful elaboration on the models of contact process can be found in |19j). 

At the other extreme, the contact rate might be independent of host density. Assum- 
ing that susceptible and infective were randomly mixed, this would lead to transmission 
following f3SI/N = (331 /{S + I). This mode of transmission is often called 'frequency- 
dependent' or 'density- independent' transmission [2n|- It is often assumed in models of 
sexually transmitted diseases because the number of sexual partners of an individual usu- 
ally depends on the mating system of the species and is weakly related to host density 

We note here that the distinction between these two modes of transmission is not 
crucially important for many problems in human diseases, because most pathogens cause 
little mortality and the total population size remains more or less constant |22j . 

Even if some demography processes are taken into account many epidemiological mod- 
els are formulated so that the infectious disease spreads in a population that has a fixed 
size N = const with balancing inflows and outflows due to birth or death or migration, 
see, e.g., |23l[22]- However, if the population growth or decrease is significant (it might be 
the case on long periods of time, e.g., for diseases transmitted vertically) or the disease 
causes enough death to influence the population size, then it is not reasonable to assume 
that the population size is constant. In this case the difference between the mass action 
transmission and the frequency-dependent transmission becomes profound p5lj23j. It is 
worth noting that the frequency-dependent transmission is considered to be more realistic 
for most human diseases |22j . 

Due to the particular form of the transmission function, (381 /{S + /), the origin 
S* = 0, / = becomes undefined in models where host extinction is a possibility, and the 
properties of the origin crucially affect the global dynamics of the system. 

As an example of such models, we rewrite a model of Chagas' disease from 
that takes the form 

-^ = {h-r- v)S + (6i(l -q) + c)I-(3 



dt S + 1 ^ 

(piq — ri — c)I + vS + [3 



dt ' " ' ' s + r 

where b, bi, r, ri, c, q, respectively, denote the birth rates of susceptible and infective 
individuals, their death rates, the cure rate and the probability of vertical transmission; 
V is the transmission rate of disease by a vector. From we have 

l(5 + /) = (6-r)5 + (6i-ri)/, 
and the total population could be either increasing or decreasing. 
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A usual way to analyze models like ([U]) is to perform the change of variables x = 
S/{S + I), y = I/{S + I). After this and using the fact that x + y = 1, it is often 
straightforward to obtain the analysis of the system in coordinates x, y, simultaneously 
keeping track of the asymptotic behavior of the total population size N [23l [25l 
However, after this transformation some of the information concerning initial variables S 
and I can be lost. For example, in variables x and y it is difficult to see that the origin 
can be simultaneously attractive and repelling, and deterministic extinction of the total 
population can be accompanied by an initial growth of the susceptible subpopulation (see 
Section m Example 3). 

Though mathematical models of pathogen transmission with function j3SI/N were 
extensively studied for years, the first case study, to our knowledge, that focuses on the 
possibility of deterministic extinction of host population was considered in [27] . A careful 
analysis of a more general model was presented in [28j. 

As an example, let us state the model of host-parasite interaction from This model 
allows for host extinction and, formulated through the basic birth and death processes, 
has the form 

J Q Q J 

a{S + J) - ail - 6)1 - cS{S + I) - [d + m)S - (3 



dt S + 1 /«\ 

dl ^ ' 

- = + a)/ -c/(S + /)+;3—. 

where a, c, d, m, a are nonnegative parameters, and 0^6'^ 1. Model ([7]) was reduced 
to a Gause-type system by means of one blow-up transformation [S, I) —>■ [u, I), u = S/I. 
We show (Section HI Example 2) that using this transformation is not enough to obtain 
a complete qualitative picture of a neighborhood of 0(0, 0). 

Systems ([6]) and ([7]) can be transformed into polynomial systems with a non-hyperbolic 
equilibrium at the origin by the time change dt ^ {S + I)dt, and the resulting systems 
of ODEs are in the class ([T]). 

It is worth mentioning that if the pathogen transmission follows mass action kinetics, 
but 5* and / denote numbers and not densities, and the total population density (here we 
speak of local spatial density) remains constant, one has to use [3SI /{S + I) expression 
to model transmission of disease 



Other models. We by no means gave a full list of models analysis of which requires 
the knowledge of the structure of VL of undefined equilibrium 0(0,0) and which can be 
reduced to ([T]). Similar models appear in many other areas of mathematical modeling in 
biology. Some mathematical models of interaction between populations of cells with a 
virus population can also be formulated in the form ([T]); e.g., one such model was applied 
to simulate anticancer therapy with oncolytic viruses [HI] (see Section HI Example 1). 
It was shown that the model possesses dynamical regimes that lead to elimination of 
cancer cells (a phenomenon known from clinical trials). Recently it was suggested that 
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immune system response is more consistent with empirical data if it is considered to be 
ratio-dependent, see, e.g., [5^ . 



3 Non-hyperbolic equilibrium of system (Hj) and struc- 
tures of its neighborhood 

3.1 Preliminaries 

Let 0(0, 0) be an isolated singular point of the vector field 

J{x,y) = P{x,y)— + Q{x,y)— (8) 
and, correspondingly, an isolated equilibrium of the system of ODEs 

(9) 

Here 

P{x, y) = P„(x, y) + P*(x, y), Q{x, y) = Qn{x, y) + Q*{x, y), 
where Pn{x,y), Qn{x,y) are homogeneous polynomials of the n-th order: 

Pn{x, y) = Pn,lX'^ + Pn-l,2a;""^2/ + . . . + P0,n+iy'', 
Qnix, y) = Qn+lflX'' + qn,lX''~^y + ■■■ + gi,n2/", 

and P*{x,y) = 0{\{x,y)\^+'), Q*{x,y) = 0{\{x,y)r+'). 

It is well known (e.g., [33]) that point 0(0,0) can be either monodromic (i.e., a focus 
or a center) or all orbits of system tending to the origin for t ^ oo or t — — oo 
(we shall call such orbits as O -orbits), approach the origin along characteristic directions. 
The latter means that any 0-orbit has definite direction 9 so that arctan(?/(t)/2;(t)) 9 
for t— »-cx)ort— >— cxD. We note that the number of characteristic directions is finite in 
the generic case. Hereafter we refer to Andronov et al. [33j for a number of results and 
notations used, but restate some of the results in the form convenient for our exposition. 

Following the results from Ch.VIII Sec. 17 of [33J we can state that any small enough 
neighborhood VL of an isolated equilibrium point of an analytic system can be partitioned 
by 0-orbits with characteristic directions into open regions, called sectors. These sectors 
can be classified into three types: parabolic sectors, hyperbolic sectors, and elliptic sectors, 
respectively; we shall call them the Brouwer sectors. The Brouwer sectors are described 
in the following figure (Fig. [T]) and their definitions are given, e.g., in [33] • Note that the 
topological equivalence of a sector to one of the sectors in Fig. [T] need not preserve the 
directions of the flow. 

More exactly, the next statement is valid. 
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(a) (b) (c) 




Figure 1: The Brouwer sectors, (a) A parabolic sector; (b) a hyperbolic sector; (c) an 
elliptic sector 

Theorem 1. There exists a neighborhood Q of a non-monodromic singular point of Ca- 
veator field, such that any neighborhood Q G Q that contains this singular point can be 
partitioned into a finite number of sectors of parabolic, hyperbolic and elliptic types. The 
number, types and the cyclic order of a position of these sectors hold with decreasing of 
the size of neighborhood Q. 

For n = 1 we have 

Pi (x, y) =piix+ po2y , Qi {x, y) = q2oX + quy . 

Equilibrium 0(0,0) of ([9]) with n = 1 is hyperbolic ii D = puqu + J9o25'20 7^ and 
Tr = pii + gii 7^ for D > 0. Point 0(0,0) has 0-orbits with characteristic directions 
if Tr^ — 4D ^ 0. If the last condition holds and additionally D ^ then 0(0,0) is a 
saddle whose neighborhood Q contains four hyperbolic sectors, or 0(0, 0) is a node whose 
neighborhood Q contains only parabolic sectors. In general, the type of the origin of 
system ([9]) in the case = 1 is determined by the calculation of eigenvalues, and there is 
a simple algorithm to infer the structure of a small neighborhood of this point. 

For n ^ 2 the situation is more complicated. Equilibrium point 0(0,0) is not hy- 
perbolic (both eigenvalues are zero). The following problem is of principal interest: to 
describe all possible topological structures of a small neighborhood Q o/O(0, 0) and asymp- 
totic behavior of orbits with variation of the system coefficients. In the present work we 
solve this problem in an important case n = 2 under some additional assumptions on (I9j). 

The problem of finding and describing the sequence of the Brouwer sectors that 
constitute fl, as well as asymptotics of 0-orbits, was studied in a number of works 
[311 ESI ESI EH EHl EH]- The main attention was paid to constructing algorithms of analy- 
sis of a non-hyperbolic equilibrium with the parameter values fixed. The general blow-up 
method, developed in [38l [391 [36] for so-called 'system with a fixed Newton diagram', 
allowed to describe topological structures of fl and compute asymptotics of 0-orbits close 
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to the singular point. Below we apply the developed methods to system with n = 2 
and present an approach to analyze the behavior of with parameter changes. 

3.2 Basic theorems 

Let us consider system ([9]) with n = 2 in neighborhood Q of an isolated equilibrium 
0(0,0), and let 

y) = xQnix, y) - yPn{x, y). 

Definition 1. We shall call vector field (IHl) {and the corresponding system of ODEs ([9])) 
non- degenerate in if 

(CI) polynomials Pn{x, y), Qn{x, y) have no common factors of the form 

ax + hy, where at least one the constants a, b is non-zero; 
(C2) polynomial F{x,y) has no factors of the form {ax + by)^, where k > 1. 

The coefficients of polynomials P2{x,y) and Q2{x,y) can be considered as the system 
parameters, and the parameter space is divided into domains of topologically equivalent 
behavior. The main results of our analysis show that each parameter domain corresponds 
to one of the following four types of phase portraits (Fig. [2]). 

Theorem 2. Let system with n = 2 be non-degenerate. The positional relationship 
of the Brouwer sectors in Vl can be of four topologically non-equivalent cases: 

(i) six hyperbolic sectors {Fig. [Ha); 

(ii) two hyperbolic sectors separated by two parabolic sectors such that one of the 
parabolic sectors is attracting (in a sense that 0-orbits tend to 0(0,0) for t ^ oo) and 
another is repelling (0-orbits tend to 0(0,0) fort —oo) {Fig. IM>); 

(iii) two hyperbolic sectors (Fig. {Mp); 

(iv) two elliptic sectors separated by two parabolic sectors such that one of the parabolic 
sectors is attracting and another is repelling {Fig. l^). 

Remarks to Theorem [21 

1. In the case considered there are always 0-orbits with characteristic directions. 

2. The boundaries between the parameter domains are given by violations of condi- 
tions {CI) and (C*2). Violation of {CI) leads to the presence of a straight line, passing 
through the origin, of non- isolated equilibrium points of ©, violation of (02) may lead 
to disappearance of characteristic directions. 

3. The topological equivalence to one of the presented in Fig. [2] cases need not preserve 
the directions of the flow. 

To describe possible asymptotics of 0-orbits of system Q we require the following 
condition: 

(03) Pn{0,y)=0^ P*{0,y) = 0, Qn{x,0) = ^ Q*{x,0) = 0. 
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Figure 2: Four possible types of the topological structure of fl of 0(0,0) for system ([9]) 
with n = 2 

Recall that curve y = f{x), /(O) = has an exponential asymptotic with positive 
power p and non-zero coefficient k if 

f{x) = kx^il + oil)), k^O,p>0. (10) 

Theorem 3. Let system ([9]) be non-degenerate. Then 

(i) An O- orbit with a characteristic direction has asymptotic f IlOp with p = 1 if and 
only if polynomial F{l,u) has a root u = k 0; if P2{0,y)Q2{x,0) ^ then system ^ 
has no O- orbits with other asymptotics; 

(ii) IfQ2{x, 0) = and (C3) holds then system ([9]) has trivial orbits y = forx > and 
X < 0; if additionally (3 = g2,i/P2,i > 1 then ([9]) has a family of exponential asymptotics 
ffTUj) with p = (3, and k ^ Q is an arbitrary constant; 

(iii) // P2(0,?/) = and (C3) holds then system has trivial orbits x = for 
y > and y < 0; if additionally {3 = Pi^2/li,2 > 1 then Q has a family of exponential 
asymptotics ( fTOl) with p = 1/13, and k ^ Q is an arbitrary constant. 

The proofs of Theorems [2] and [3] are given in Appendix. 
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3.3 An algorithmic approach to the analysis of the structure of 
the origin of system ([TI) 

The results from the previous section show what can be expected in neighborhood Q of the 
origin of system ([T]). Here we summarize a practical recipe to construct a phase-parameter 
portrait of ([I]), i.e., we show how to analyze a particular system. The presented algorithm 
follows directly from the proof of Theorem [2] (see Appendix). 

Due to the homogeneous form of system ([1]) the analysis of an isolated equilibrium 
0(0, 0) can be formulated in terms of functions depending on one variable. Let F{x, y) = 
xQ2{x, y) — yP2{x, y). We need to consider two pairs of polynomials 

(i) P2(l,n), Fi(m)=F(1,m), 

and 

(ii) Q2{vA), F2{v) = -F{v,l). 

The conditions of non- degeneracy (CI) and (C2) take the form: 

(CI') The pair of polynomials P2(l,^), ^2(1,^), as well as the pair of polynomials 
P2(^5 1)5 <52(^^5 1)5 have no common roots; 

(C2') The polynomials Fi{u) and F2{v) have no multiple roots. 

The analysis relies on finding the roots of Fi{u) and F2{v). There can be five different 
cases, where we list only the roots that are necessary for the subsequent calculations: 

(i) Fi{u) has three real roots Ui, i = 1,2, 3; 

(ii) Fi{u) has two real roots Ui 0, i = 1,2, and zero is a root of F2{v); 

(iii) Fi{u) has two real roots Ui = 0, U2 7^ 0, and zero is a root of F2{v); 

(iv) Fi{u) has one real root; 

(v) Fi{u) has no real roots, and zero is a root of F2{v). 

The lines y = UiX and, if zero is a root of F2{y), x = divide VL of 0(0,0) into six 
(cases (i)-(iii)) or two (cases (iv)-(v)) sectors with six or two branches. These sectors can 
be of a hyperbolic, elliptic, or parabolic type. Each line consists of two branches divided 
by the point 0(0,0). 

We introduce the following notations. For the branches of lines y = UiX we consider 

for the branches of line x = v = we consider 

Xliv)=Q2iv,l), K{v) = F^{v). 

Note that we need to use functions F2{v) and Q2{v, 1) only in cases (ii), (iii) and (v). 
We shall call numbers A", Ag, A" and Ag the branch characteristics. 
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Figure 3: Placing of the Brouwer sectors depending on the type of the sector branches {h 
- hyperbohc branch, p - parabohc branch) 

Definition 2. We shall call a branch of a sector in the phase space {x,y) of system ([T]) 
hyperbolic if the inequality X1X2 < (or XIX2 < 0) holds, and parabolic z/ A^Ag > (or 
X1X2 > 0) holds, where A", Aj (or A^', Ag) the branch characteristics. 

Let be a sector in the state space of system ([T]) with a vertex at 0(0, 0) composed 
by branches of the hnes y = uix and y = U2X (or one of the hnes can be a; = 0), and 
Vq = V f]Q. The following proposition, proved in Appendix, holds. 

Proposition 1. contains 

(i) a hyperbolic sector if both of its branches are hyperbolic [Fig. 0a); 

(ii) an elliptic sector if both of its branches are parabolic {Fig. H)); 

(iii) a parabolic sector {or its part) if one of the branches is hyperbolic and another is 
parabolic {Fig. \Ep)- 

The direction of the phase flow on the sector branches is determined by the sign of 
Ai. If Xi{ui) < then the phase flow goes to 0(0, 0) for a; > and from 0(0, 0) for a; < 
(see, e.g., the line y = Uix in Fig. [3^); if A"('Ui) > then the phase flow goes from 0(0, 0) 
for X > and to 0(0, 0) for a; < (the line y = U2X in Fig. EK)- If A^(0) < then the 
phase flow goes to 0(0,0) for ?/ > and from 0(0,0) for y < 0; if X\{Q) > then the 
phase flow goes from 0(0, 0) for ?/ > and to 0(0, 0) for y < 0. 

The conditions that determine the boundaries of topologically non-equivalent domains 
in the parameter space now can be reformulated in terms of the branch characteristics. 
Namely, 
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4 Examples 



In this section we give a number of examples of analysis of the complicated non-analytical 
equilibrium point in biological models. Mainly we restrict our attention to M^, but em- 
phasize that in some cases it is necessary to analyze not only the behavior of the state 
variables in but also the behavior in adjacent areas. 

Example 1. {Mathematical model of anticancer treatment with oncolytic viruses) On- 
colytic viruses that specifically target tumor cells are promising new therapeutic agents 
[40j. The interaction between an oncolytic virus and tumor cells is highly complex and 
nonlinear. Hence, to precisely define the conditions that are required for successful ther- 
apy by this approach, mathematical models are needed. Our model [31] was formulated 
through the incorporation of frequency-dependent mode of virus transmission into the 
model of Wodarz [H]. The model, which considers two types of tumor cells growing in 
logistic fashion, has the following form (non-dimensional variables and parameters are 
used) : 



^ = x(l - (x + y)) - /3-^ , 
at X + y 

— = 7^(1 -{x + y)) + 13 — ■ by. 

dt X + y 



fill 



X + y 

Here x, y are (scaled) sizes of non-infected and infected tumor cell populations, respec- 
tively; 7, /5, 5 > are non-dimensional parameters (/3 takes stock of the transmission 
rate, and 5 describes the virus cytotoxicity). 

After the time change dt ^ {x + y)dt we obtain the system of ODEs in the form (Q. 
We have 

P2{x, y)=x^ + {l- P)xy, Q2{x, y) = (7 - 5 + p)xy + (7 - S)y\ 
The polynomial 

Fi{u) =Q2{l,u)-uP2{l,u) = (7-5 + /3-l)M(n + l), 

if 7 — 5 + /? 7^ 1, has two roots Ui = and U2 = —1. Polynomial F2{v) also has two roots, 
one of which ■0 = (another one is ^2 = 1/^2 and we do not need to consider it). is 
confined by the branches of the lines y = Uix = and x = v = 0. This means that the 
root U2 is redundant for our analysis as far as we are concerned with the system behavior 
in M^. The branch characteristics are 

A^(0) = 1, A^(0) = 7-5 + /3-l, 

and 

A^(0)=7-5, A^(0) = -A2"(0). 
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Figure 4: Phase-parameter portrait of a neighborhood of the origin of system ( fTTl) given 
as a cut of positive parameter space for an arbitrary fixed value of 7 > 0. The bifurcation 
boundaries are ai = {(5,7) : S = 7}, and a2 = {(5,7) : P = S + 1 — •y} {h - hyperbohc 
branch, p - parabohc branch) 

Analyzing the sign of these expressions allows us to completely describe a small neigh- 
borhood of 0(0, 0) (Fig. Hland the following theorem). 

Theorem 4. For different positive values of parameters 5, (3, and 7 there exist three types 
of topologically different structures of the neighborhood fl+ = Qf] of point 0(0, 0) {and, 
correspondingly, three topologically different phase portraits of system ( JTTl) ).- 

(i) a repelling parabolic sector {domain I in Fig. [7]) for the parameter values 5 < 7. 
The phase curves of the system that tend to 0(0, 0) are of the form 

y = Cx'+'-^{l + o{l)), (12) 

ifP>6 + l — J, and 

y = Cx^^-^^/^^-^\l + o{l)), (13) 

z/ /5 < (5 + 1 — 7; here C ^ is an arbitrary constant; 

(ii) an elliptic sector {domain II in Fig. ^ composed by trajectories tending to 0(0, 0) 
as t ^ 00 {with asymptotics given by (fT3l)). as well as with t ^ —00 {with asymptotics 
given by ( fT2l) ) if 6 > •y and /5 > 5 + 1 — 7; 

(iii) a saddle sector {domain III in Fig. ^ for the parameter values 5 > 7 and 
< 5 + 1 -7. 

Remark to Theorem [41. Note that only one of two possible arrangements of the 
sector branches is shown in domain / of Fig. HI It is possible that the branch x = 0, y > 
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may be parabolic, and the branch y = 0, x > may be hyperbohc. In both cases, though, 
we obtain that 0(0, 0) is a repelhng node, and fi+ contains a parabohc sector. 

The most beneficial domain of parameter values is domain II in Fig. H] (an elliptic sec- 
tor). This domain corresponds to the total elimination of both cell populations (infected 
and uninfected) regardless of the initial conditions (this dynamical regime is observed in 
experimental studies, e.g., |l2])- However, on its way to extinction, the overall tumor size 
X + y can reach rather high values (which, with the parameters fixed, crucially depend 
on the initial conditions). This indicates that we must not only identify the conditions 
that favor tumor elimination, but also develop the optimal strategy to infect the initial 
tumor. In the framework of the considered model this can be done using the information 
on asymptotics ( fT2|) and ( fT3|) (see [3lj). 

Example 2. {Parasite-host interaction model) Hwang and Kuang [21] formulated a 
host-parasite model allowing for host extinction. The model, which is a non-dimensional 
version of (CI), is of the form 

— = x + y— [l— d)y — x{x + y) — ox — s 



dt x + y 

dy , . , xy 

_- = -(5 + r)y - y{x + y) + s — — , 
dt x + y 

where all parameters are non- negative, and ^ ^ 1. Solutions of fll4p are considered 
in R^. Using the transformation v = x/y, Hwang and Kuang reduced model (JT^ to a 
Gause-type system which was completely studied. 
After the time change — > (x + y)dt we obtain 

P2(x, y) = il- S)x^ + {l-6-s + e)xy + 
Q2{x, y) = {s-6- r)xy -{5 + r)y^, 

F^{u) = {-r + s - l)u+ {-r - 1 + s - 9)u'^ - 9 u^. 

The polynomial Fi{u) has three roots 

Ml = 0, M2 = -1, Mg = (S - r - l)/6', 

if Us 7^ 0, —1 (case (i) in Section 

We note that if 5 > 1 then, from (fT^ . ^{x + y) < follows. Thus, we shall consider 
only the case when < 6 < 1. 

The state space may be contained in two sectors, so in this case we need to calculate 
the branch characteristics for the branches of three lines y = UiX. We obtain 

Xliu,) = 1-5, A^(mi) = s - r - 1, Xi{u^) = s, A^(n2) = -{s-r-l + 0), 
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Figure 5: Phase-parameter portrait of a neighborhood of the origin of system (JT^ given 
as a cut of positive parameter space for arbitrary fixed values of < < 1 and < 5 < 1. 
The bifurcation boundaries are ai : A"('U3) = and a2 '■ A2 = (/i - hyperbohc branch, 
p - parabohc branch) 

and 

Xi{u^) = ~{5 + r)A + s, A^(m3) = -(s - r - l)A, 

where A= [s — r — 1) / 9 + 1. 

Using A2('Ui) = and \i{u^) = as bifurcation curves we can draw the complete 
parameter portrait in Vt^ = Qf]M.'^ of the origin (Fig. [5]). 

In general, we obtain 

Theorem 5. Let < 5 < 1. For different positive values of the parameters there exist 
three types of topologically different structures of the neighborhood i7+ of point 0(0,0) of 
system iHM : 

(i) An elliptic sector and a part of an attracting parabolic sector for the parameter 
values Xiiu^) > {domain I in FigurelS^. The phase curves of the system that tend to 
0(0, 0) when t —00 are of the form 

y = Cx(^-'-5)/(i-5) (1 + 0(1)); (15) 

(ii) A part of a hyperbolic sector and a repelling parabolic sector for the parameter 
values Xiius) < and A2('Ui) > {domain II in Figurel^. The phase curves of the 
system that tend to 0(0, 0) when t —00 are of the form 

(iii) A part of a hyperbolic sector for the parameter values X2{ui) < {domain III in 
Figure 0) . 
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Remark to Theorem O. In Fig. Elthe line X'^i''^'^) = is also shown (0:3). Analysis of 
mutual placing of the branches of the Brouwer sectors shows that the topological structure 
of does not change when we cross 03. 

Thus we have proved the existence of an elliptic sector in model fll4p . which was 
overlooked in the original analysis of Hwang and Kuang [22] ■ This dynamical regime 
corresponds to the initial growth of the total population size, and the eventual population 
extinction can be preceded by a relatively normal population evolution. 

Example 3. [Model of Chagas' disease [23l [26] ) The model of Chagas' disease from 
[231 Eg has the form 

^ = (b — r — v)x + (61(1 — q) + c)y — P ""^^ 



dt x + y 

dy \ n 

— = [hiq - ri - c)y + vx + [3 — ■ — , 

dt x + y 

where all the parameters are nonnegative and ^ g ^ 1 (see also Section [21 system 
After the time change dt ^ {x + y)dt we obtain a model in the form ([T]), where 

P2{x, y) = (b — r — v)x'^ + (6 + 61 — r — biq — v — k)xy + (61 — biq)y'^, 
Q2{x, y) = vx^ + {biq -ri + v + k)xy + {biq - ri)y'^. 

In the following we assume that 6 — r > and a = 6 — 61 + ri — r > 0. 
The polynomial 

Fiiu) = {u + l)(-&i(l - q)u^ + (A; + w - a - 61(1 - q))u + v) 

has three roots 

, . A±JA^ + Avbi{l-q) 

= = mi^) ' 

where A = k + v — a — bi{l — q) . We have U- < and u+ > ii v > 0. As in the previous 
example, may be contained in two sectors, so we need the branch characteristics of 
all the branches: 

A"(mi) = k, X2{ui) = a - k, 

and 

Ai(m±) = (6i-ri)(M± + l) + a. 

We do not explicitly present A2(m±), but it is easy to see the signs of these expressions. 
We know that Fi{u) is a polynomial of the third order, Fi{oo) = —00, and Fi{u) has 
one positive and two negative roots. This implies that Af («+) = F{{u+) < for any 
parameter values. The sign of X2{u-) can be either negative or positive depending on the 
sign of A2('Ui) (they have the opposite signs). 

Analyzing the branch characteristics we obtain the following parameter dependent 
structure of = ^ fl ^+ '^^ origin. 
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(a) b^-r^>0 
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Figure 6: Phase-parameter portrait of a neighborhood of the origin of system (1161) given 
as a cut of positive parameter space for (a) &i — ri > and (6) 6i — ri < 0. The bifurcation 
boundary is ai = {(6, 6i, r, ri, k, v, q) : k = — ^ + {h - hyperbohc branch, p 

- parabohc branch) 

Theorem 6. For different positive values of the parameters there exist three types of 
topologically different structures of the neighborhood Qj^ of the origin of system (11 6 1) ; 

(i) A part of a repelling parabolic sector for the parameter values 6i — ri > and 
Xi{u ) < {domain I in FigurelBi); 

(ii) A part of a hyperbolic sector and a part of a repelling parabolic sector for the 
parameter values 6i — ri > and \i{u^) > or bi — ri < and A"(m_|_) > {domain II 
in FigurelB^; 

(iii) A part of an elliptic sector and a part of an attracting parabolic sector for the 
parameter values 6i — ri < and A"(u+) < {domain III in FigurelB^. 

Theorem [6] gives some additional information about the behavior of the state variables 
in comparison with the original analysis [23]. Most importantly we have found that pop- 
ulation extinction occurs when neighborhood Q of 0(0, 0) contains an elliptic sector, only 
part of it belonging to M^. Depending on the parameter values, though, this elliptic sector 
can be almost entirely in (domain 1 1 lb in Fig. [U]). In other words, the population 
extinction can be essentially non-monotonous, and the total population size can reach 
relatively high numbers before vanishing. 

The bifurcation curve cti corresponds to the presence of a line of non-isolated equi- 
librium points of the form {x*,y* = —{b — r)/{bi — ri)x*). In the case 6i — ri < this 
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line belongs to M^, and the system behavior significantly changes when ai is crossed. 
If parameter values are in Domain He in Fig. the population size goes to infinity; if 
parameter values are in Domain /// the population becomes extinct with time; while if 
we are on ai, the line of non-isolated equilibrium points is an attractor, the asymptotic 
size of the population is finite and non-zero, and this asymptotic size depends on the 
initial conditions. 

It is worth mentioning that our local analysis gives the full description of the global 
behavior of system (fT6l) in a finite part of the phase plane. 

Due to Theorem [2] we can add any terms of order 2 or higher to the right side of (fT6ll , 
and the parameter portrait of f2+ will not change. For example, it is reasonable to assume 
a logistic regulation of the population p3]. 

5 Conclusions 

In this paper we presented qualitative methods to analyze a wide class of models of bio- 
logical populations and communities. These models are characterized by the presence of 
complicated singular equilibria. A significant number of predator-prey, host-parasite, epi- 
demiological, etc., models, which were recently suggested in the literature, belong to this 
class. The main mathematical peculiarity of these models is their non-analytic structure 
of the origin, which is a result of modelling attempts to better reflect characteristics of 
the simulated processes at small population sizes. 

Removable singularity peculiar to these models turns into an analytic equilibrium by a 
change of independent variable, the models become well-defined at this equilibrium. How- 
ever, the linearization approach fails to infer topological structure of a small neighborhood 
of this point, the origin becomes a non-hyperbolic point, i.e., it has zero eigenvalues for 
any parameter values. 

Methods of analysis of non-hyperbolic points for two-dimensional ODE systems with 
fixed values of the system coefficients were developed earlier [SU EHl EH] • A neighborhood 
Q of the origin of a smooth second order system is divided into a finite number of sectors 
with different phase behaviors: parabolic, hyperbolic, and elliptic ones. The parabolic 
and hyperbolic sectors can be realized in a neighborhood of a hyperbolic equilibrium, 
while the elliptic sector is intrinsic only to a non-hyperbolic equilibrium point. 

We presented a complete description of the structure of fl and gave asymptotics of 
0-orbits for generic systems whose Taylor series at the origin start with second order 
terms. An exact algorithm is suggested to analyze the structure of the origin when 
the system parameters vary. The theory and methods are illustrated by applying the 
developed algorithm to the phase-parameter analysis of some existing models that fall 
into the class of ODEs ([T]). In particular, a model of anticancer treatment with oncolytic 
viruses, suggested in [31], is analyzed in a neighborhood of the origin, the parameter 
domain that corresponds to the tumor elimination is identified; a mathematical model of 
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parasite-host interaction from [29j shows that the parasite population can drive the host 
to extinction with a preceding population outbreak; a model of Chagas' disease [23] also 
demonstrates the effect of essentially non-monotonous population extinction. 

The main attention is paid to the problem of existence and practical finding of elliptic 
sectors in the phase plane. An elliptic sector represents a family of homoclinics to equilib- 
rium, which means that orbits, starting in this sector, tend to the origin both for t — >■ cx) 
and t — s> — oo. Although such type of dynamic behavior was known in the theory of 
dynamical systems for a long time, only relatively recently it was discovered in biological 
models. Interpretations of this behavior may be very fruitful. It can be considered as 
a specific form of community extinction accompanied by preceding population outbreaks 



From the perspective of biological conservation, when population coexistence is desir- 
able, the parameter domains that are characterized by the presence of an elliptic sector 
should be avoided and their boundaries have to be considered 'dangerous', since both pop- 
ulations go to a rout to extinction after crossing them. From the perspective of biological 
control, when population extinction of one or more populations is desirable, this param- 
eter area is the most interesting. In this case both populations are driven to extinction 
deterministically, an outcome that prey-dependent models are unable to produce. Let us 
emphasize that such peculiarity of the system behavior may be of vital importance. For 
example, it is the case in modeling anticancer therapy with oncolytic viruses. This regime 
demonstrates a possibility of complete eradication of the tumor cells [31j. 

Similar regimes may exist in more complex and realistic models with dimensions ex- 
ceeding 2. Examples of such systems include epidemiological models with the number of 
subpopulations more then two (e.g., [26l l22]). and population models with three or more 
trophic levels (e.g., [i3l HI] ). 

It is our hope that the methods and algorithm presented here could help not to over- 
look, as it sometimes happens, the important regime of deterministic extinction given by 
the presence of an elliptic sector. 

6 Appendix 

We consider the vector field 



(see |I31[28]). 




(17) 



and the corresponding system of differential equations 





where 



P(x, y) = Pn{x, y) + P*(x, y), Q{x, y) = Qn{x, y) + Q*{x, y). 
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Here Pn{x,y), Qn{x,y) are homogeneous polynomials of the n-th order, and 

P*{x,y)=0{\{x,y)r'), Q*{x,y) = 0{\{x,y)r'). 

We assume that the origin is an isolated singular point of ffTTj) : P(0, 0) = Q{0, 0) = 0, 
and let F{x,y) = xQn{x,y) — yPn{x,y). We also assume that (fT7|) is non-degenerate, i.e, 
satisfies 

(CI) polynomials Pn{x,y), Qn{x,y) have no common factors of the form 

ax + by, where at least one the constants a, b is non-zero; 
(C2) polynomial F{x, y) has no factors of the form {ax + by)'^, where k > 1. 

To analyze vector field ffTTj) in a small neighborhood Q of the origin we apply the 
blowing-up transformations (x, y) {x, u) 

X = X, u = y/x, X ^ 0, (19) 

and {x,y) {v,y) 

y = y, v = x/y, y ^ 0. (20) 
These transformations have the following properties (Ch. 9, Sec. 21 of [33]). 

Proposition 2. 

(i) Transformation ( |T9i) defines a topological mapping of the x-slit {x,y) -plane onto 
{x,u)-plane. Points of the first (second, third, fourth) quadrant in the {x,y)-plane are 
mapped respectively onto points of the first (third, second, fourth) quadrant in the {x,u)- 
plane. The inverse transformation is defined on the axis x = of the {x,u) -plane, but 
maps it onto a single point (0,0) of the {x,y) -plane. 

(ii) Transformation ( l20l) defines a topological mapping of the y-slit {x,y) -plane onto 
{v,y) -plane. Points of the first (second, third, fourth) quadrant in the {x,y) -plane are 
mapped respectively onto points of the first (second, fourth, third) quadrant in the {v,y)- 
plane. The inverse transformation is defined on the axis y = of the {v,y) -plane, but 
maps it onto a single point (0,0) of the {x,y) -plane. 

First we apply transformation ffTIJl) . After this transformation and the time change 

dt x^'-^dt (21) 



we obtain the system 



— = xPnil, U) + Gi{x, U), 

du (22) 
= Fi{u) + G2{x,u), 
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where 



Fi{u) = u), Gi{x, u) = P*{x, ux)/x^ \ G2{x, u) = {Q*{x, ux) — P*{x, ux)u)/x'^ . 

Let M be a root of Fi{u). Then (0, u) is an equihbrium point of ([22]) with eigenvalues 
P„(l,u) and F[{u). If vector field ffTTj) is non-degenerate this equilibrium point is hyper- 
bolic, i.e., Pn{l,'u)F[{u) 7^ (hereafter the prime denotes the derivative). Summarizing, 
we obtain 

Proposition 3. Equilibrium point (0,u) of ([22]) is a saddle if Pn{^,u)F[{u) < and a 
node if Pn{l,u)Fl{u) > 0; this node is a sink if Pn{l,u) < 0, F[{u) < and a source if 
P„(1,m) >0, F[{u)>Q. 

According to Proposition [2] the point (x = 0, ?/ = 0) is transformed to w-axis in (x, u)- 
plane; w-axis consists of the orbits {x = 0, < m < Uj+i} and of equilibria (0,'Uj), where 
Ui are the roots of Fi{u). 

Knowledge of the order of the equilibrium points of ( l22i) on w-axis allows us to infer 
the structure of the Brouwer sectors in a small neighborhood of 0(0,0) in (x, ?/)-plane 
(except, perhaps, close to the line x = 0). To specify possible characteristic directions of 
0-orbits we need the following theorem (compare to Th. 64 from [23] )• 

Theorem 7. 

(i) Any 0-orbit of system ( |T8i) is either a spiral or tends to 0(0,0) with a definite 
tangent, i.e., has a characteristic direction; 

(ii) // at least one 0-orbit is a spiral then all orbits in some neighborhood o/ 0(0,0) 
are spirals; 

(iii) If polynomial F{x,y) = xQn{x,y) — yPn{x,y) ^ identically, then coefficients k 
of all tangent lines y = kx to 0-orbits are the roots of the polynomial Fi{u) = F{l,u) 
{except, perhaps, the tangent line x = corresponding to the root v = of polynomial 
F2{v) = -F{vA)). 

The structure of non-degenerate vector field ( JT71) between the lines, corresponding to 
neighboring roots of Fi{u), is described by the following proposition. 

Proposition 4. Let Ui < Ui^i be neighboring roots of Fi{u), i.e., there are no other roots 
of Fi{u) between them. 

(i) If points {0,Ui), (0,iij+i) are saddles, then neighborhoods, of the origin of (|T8l) has 
two hyperbolic sectors, one for x > and another for x < 0, whose separatrixes are the 
curves 

af. y = Uix{l + o(x)), a^+i : y = Mi+ix(l + o(l)); 

(ii) If points (0,-Uj), (0,'Uj_|_i) are nodes, then neighborhood Q of the origin of ([TH]) has 
two elliptic sectors, one for x > and another for x < 0. These sectors contain 0-orbits 
whose asymptotics are ai and aj+i; 
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(iii) // point (0, Ui) is a saddle and point (0, Mj+i) is a node, then neighborhood Q of 
the origin of ffT^ has two parabolic sectors, one for x > and another for x < 0. These 
sectors contain O -orbits whose asymptotic is Oj+i. 

Proof. The points (0, Uj) and (0,{ij+i) are equilibria of (1221) with eigenvalues 
and 

From the form of the Jacobian matrix it follows that u = Ui and u = Mj+i are the 
characteristic directions of orbits tending to (0,-Uj) and (Oj-Uj+i) respectively, i.e., these 
orbits have asymptotics 

7i: M = Ui{l + o(l)), 7i+i: u = Ui+i{l + o(l)). 

(i) Let both equilibria be saddles (see, e.g.. Fig. O i = 1). We have X\X\ < 0; 
Ai+U^+^ < 0, and, due to (C2), A^A^+^ < 0. It implies that AiAi+^ < 0, i.e., the flow along 
curves 7^ and 7j+i has opposite directions (at least for small x). Let point {xo,Uo), where 
Ui < uq < Ui+i and xq is small enough, belong to orbit Lq of ( l22l) . Then Lq approaches 
7j+i when t increases and 7^ when t decreases (Fig. ^ i = 1); Lq forms a hyperbolic-like 
curve with asymptotics 7, and 7i+i. Due to continuity any orbit passing through a point 
(x, m) close to (xo,|/o) is a similar curve. Applying Proposition [5] we obtain a hyperbolic 
sector in coordinates x, y. 

(ii) Let both equilibria be nodes (see Fig. [HI ^ = 2). Similarly to the above the flow 
along curves 7i and 7j+i has opposite directions, one node is a sink and another one is 
a source. According to Andronov et al. Ch. 4, Th. 20 [33] infinite number of orbits 
have asymptotics 7^ and 7i+i. Taking orbit Lq as above we obtain that Lq tends to the 
source for t ^ 00 and to the sink for t —00 forming a parabolic-like curve. The curves 
containing points {x,y) close to {xo,yo) have the same form. Applying Proposition [2] we 
obtain an elliptic sector in coordinates x,y. 

Case (iii) is considered similarly to the previous ones. □ 

Remark to Proposition HI Here it is worth noting that we also need to specify the 
direction of the vector field ffTTl) with respect to the origin. When we speak of the direction 
we mean that the flow either goes towards the origin or from the origin. For small positive 
X the direction of vector field ( |T7I) on the line with asymptotic y = ux{l + o{l)) coincides 
with the direction of the vector field corresponding to (!22|) on m = u{l + o(l)) (i.e., either 
both flows go towards corresponding equilibria or from them); for small negative x these 
directions are the same if n is odd and opposite if n is even. This simple fact follows 
from the time change (12T1) and should not be forgotten when analyzing the direction of 
the vector field. 
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The number of equilibria of ([221) of form (0, u) is equal or less than n + 1. Denote 
Ug and rin the numbers of saddles and nodes respectively. 

Proposition 5. 

Us ^ n + 1, Un < n + 1. 

Proof. Let -Uj+i be neighboring roots of Fi{u) and A*]^ = P„(l,'Uj), = F[{ui) are the 
eigenvalues of equilibrium (Oj-Uj). Due to (C2) we have that A2A2'''^ < 0. The polynomial 
P„(1,m) can change sign at most n times. Polynomials Fi{u) and P„(1,m) have opposite 
signs for large u. We have to show that the sequence {AJA2, . . . , A"'''^A2'''^} cannot have 
all the elements positive. If, e.g., -Fi(oo) = +00 then F[{un+i) > 0, where Un+i is the 
largest root of Fi{u). If we suppose that Pn{l,Un+i) > 0, it means that Pn{l,u) has 
already changed its sign at least once and we have n — 2 possible sign changes for the rest 
n roots. Here is a contradiction, and we cannot have n„ = n + 1. Similar arguing shows 
that the sequence {AJA2, . . . , XI^^X^^^} can have all the elements negative. □ 



Blowing-up transformation (1191) allows us to explore the behavior of vector field (|Tj 
anywhere except for the axis x = 0. The behavior of (|T71) close to y-axis can be investi- 
gated with the help of the second blowing up transformation (|20l) . 

After the time change 

dt y^'-^dt (23) 



we obtain the system 



dv ^ , . , 
= F2{v)+H,{v,y), 



dt 
dy 

dt 



(24) 

yQn{v, 1) + H2{v,y), 



where 

F,{v) = -F{v,l), H^{v,y) = {P*{vy,y) - Q*{vy,y)v)/y^, H,{v,y) = Q* {vy , y) / y^~\ 

According to Proposition [2] systems (|T8l) and (1241) are equivalent in all the points y 0; 
the point {x = 0,y = 0) in transformed into axis v in (f , ?/)-plane; the axis v consists of 
the orbits {vi < v < -Oj+i, y = 0} and of equilibria {vi, 0), where Vi are the roots of F2{v). 

There is a strict correspondence between the number and types of equilibrium points 
of ([22]) on M-axis, and the number and types of equilibrium points of ([21]) on f-axis. 



Proposition 6. Let vector field ( JTTl) be non-degenerate, polynomial Fi{u) have a root 
u ^0, and A^, Ag be the eigenvalues of the singular point (0,-u) of system ( I22l) . Then 

(i) polynomial F2{y) has a root v = 1/u; 

(ii) the eigenvalues X\, A2 of the singular point {v, 0) of system ([21]) satisfy the relations 
A^ = AV^^'S A^ = A^/m"-^ 
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Proof, (i) We have Fi{u) = F{l,u), F2{v) = -F{v, 1) and 
F{x,ux) = x'^'^'^ Fi{u) , F{vy,y) = —y 




Fi{u) = implies F{x,xu) = for any x. Consider 

F{vy,uvy) = 1 



which means that I/m is a root of F2{v). 
(ii) We use the following notations: 



We have = -nP„(l,u) +Q„(l,u). Thus, A^' = P„(l,u) = Qn{l,u)/u = 

u'^Qnil/u, 1)/U = M"-lg„(l/M, 1) = U^~'X\. 



Next, we can write ^(l,^) = u"'^^F{l/u, 1) for any u ^ 0. Taking the derivatives of 
this relation and putting u = u we obtain 



Corollary 1. The topological type of equilibrium point {0,u), u ^ 0, of system fl22l) 
coincides with the topological type of equilibrium point 0) of system ^2M . that is, 
both points are saddles or nodes simultaneously. 

The only root of F2{v) that cannot be described with the knowledge of types of the 
roots of Fi{u) is, due to Proposition [6l v = 0. 

Proposition 7. Let system ffTSl) satisfy (CI) and (C2) and have an isolated singular 
point at the origin. Then the point [v = 0,y = 0) is an equilibrium of fl2^ if and only if 
F2{v) has a rootv = 0. It is a saddle if F2{0)QniO, 1) < and a node if F2{0)QniO, 1) > 0; 
this node is a sink 2/^2(0) < 0, Qn{0, 1) < and a source 2/^2(0) > 0, Qn{0, 1) > 0. 

Obviously, an analogue of Proposition H] holds for blowing-up transformation ( l20l) . 
The only thing that should be emphasized here is the simultaneous directions of vector 
field ffT7|) and the vector field corresponding to ([21]) (see the remark to Proposition Hj). 

The lines y = UiX and x = Viy, where Ui and Vi are the roots of polynomials Fi[u) 
and F2{v) respectively, divide a neighborhood of the origin of ( |T8l) into sectors, whose 
boundaries are the branches of these lines (each line consists of two branches). 

Below in the proof of Proposition [1] we use the notations introduced in Proposition [61 

We shall call the branches, corresponding to the line y = UiX, of a sector in phase 
space (x, y) hyperbolic if for the eigenvalues of the equilibrium (0, Ui) the inequality 
Xi{ui)X2{ui) < holds, and parabolic if X'^{ui)X2{ui) > holds. An analogous defini- 
tion applies to the branches corresponding to the line x = Viy. 



X^2 = i^:(l, u) = {n + 1)^"F(1/^, 1) - ^"-^F:(t;, 1) 



v=l/u 



□ 
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Proof of Proposition [H If a sector belongs to the half plane x > (this also implies that 
there is a complementary sector for x < 0) or ?/ > 0, the assertion follows from Proposition 
Hlor from its analogue for the blowing-up transformation ( !20l) . There are two cases which 
are not covered by Proposition HI 

(i) A sector is not contained in half plane x > 0, as well as in ?/ > (e.g., all the roots 
of Fi{u) satisfy the condition -Uj > 0); 

(ii) A sector coincides with the first (or second) quadrant of the phase plane, i.e., u = 
and '0 = are the roots of polynomials Fi{u) and F2{v). 

We start with the case (i). 

Due to the assumptions we have that x = is not a characteristic direction and the 
sector under consideration is confined by the branches of lines corresponding to the least 
and the largest roots of Fi{u), which we denote u and U. 

The main idea of the proof is exactly the same as the one used in Proposition HI We 
need to show that the vector field is coordinated. The required statement then follows 
from the continuity of the vector field. 

Let n be even. The polynomial Fi{u) can have 2k + 1 roots, which implies that 
X2{u)X2{U) > 0. In the proof of Proposition H] we had the opposite inequality for two 
neighboring roots, but considered the sectors, both branches of which lay in the half plane 
X > 0. Due to the time change (1211) for even n the direction of the vector field of system 
(|T8l) on the asymptotic y = ux is opposite to the direction of the vector field of system 
(1221) on the asymptotic u = u for x < 0. Repeating the proof of Proposition H] and taking 
into account the direction of the vector fields (e.g., if both roots correspond to saddles of 
system this yields that A"(ii)A"(t/) > and the vector field on the branches of the 
sector has opposite directions in respect to the origin) we prove the claim. 

The case of only one root of Fi (u) is easily included in the above reasoning. 

If n is odd, polynomial Fi{u) has 2k roots. We assume here that k > (see also 
general remarks at the end of Appendix). It implies that A2('u)A2(t/) < 0. Now the 
direction of the vector field on the branch ux coincides with the direction of the vector 
field on the line u = u of system (|22l) . and using the reasoning from Proposition H] we 
prove the statement. 

(ii) We consider the case when a sector coincides with the first quadrant. The axes 
are asymptotics to the orbits starting in the first quadrant because there are no other 
asymptotics of the form y = UjX, where Ui > 0. From Proposition [2] it follows that 
systems ([22]) and ([21]) are equivalent inside the first quadrant and we can consider them 
simultaneously. 

We have that m = is a root of Fi{u) and t) = is a root of F2{v). This and non- 
degeneracity of the vector field (fT71) mean that Taylor series of Fi{u) and F2{v) start from 
the terms of order n. 

If n is even, then Fi{u) can have 2k roots, one of them is zero. If we denote the least 
root of Fi{u) as u, we obtain A2(0)A2('?i) < 0. The root v = l/it is the largest root of 
F2{v) (except for zero), and, from Proposition [6l we have A2(u)A2('0) < 0. This yields 
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Figure 7: Correspondence between the phase plane (x, u) of system (l22i) and neighborhood 
Q of the origin of system (fT8|) . The case of six hyperbohc sectors 

that A^(0)A^(0) < 0. 

If n is odd, then again similar arguing shows that A2(0)A2(0) < 0. The rest of the 
proof repeats the proof of Proposition HI □ 

Now we are ready to prove the basic theorem. 

Proof of Theoreml^ We have n = 2 and polynomial Fi{u) can have three, two, one or no 
real roots. 

First we consider the case of three real roots of Fi{u). Note that we can include 
the case when one of the roots is zero (in this case, due to Proposition [HI F2{v) has 
two non-zero roots). It is possible, generally speaking, seven cases of mutual order of 
the equilibrium points of ( l22l) on w-axis: (1) three saddles; (2) node-node-saddle; (3) 
node-saddle-node; (4) saddle-node-node; (5) saddle- saddle-node; (6) saddle-node-saddle; 
(7) node-saddle-saddle. The case when there are three nodes cannot be realized due to 
Proposition O 

Placing the Brouwer sectors between the neighboring roots is accomplished with the 
help of Proposition [H 

Thus we have that case (1) corresponds to 6 hyperbolic branches and 6 hyperbolic 
sectors (Fig. [7]); cases (2)-(4) correspond to two hyperbolic and four parabolic branches 
or, respectively, two elliptic sectors and two parabolic sectors (Fig. [H]), here four parabolic 
sectors described by Proposition [T] merge into two; cases (5)-(7) correspond to four hy- 
perbolic and two parabolic branches, or to two hyperbolic and two parabolic sectors (Fig. 

ED. 

Now let Fi{u) have only one real root. We can include the case u = due to (CI) and 
(C2) {F2{v) has no real roots in this case). This root can correspond to a saddle or node 
equilibrium point (0,^) of (1221) . The first case agrees (Proposition [T]) with the presence of 
two hyperbolic sectors (Fig. [10]); the second one agrees with the presence of two elliptic 
sectors (Fig. [TTi) . 
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Figure 8: Correspondence between the phase plane (x, u) of system ( l22i) and neighborhood 
Q of the origin of system fll8p . The case of two parabohc and two elhptic sectors 




Figure 9: Correspondence between the phase plane (x, u) of system and neighborhood 
Q of the origin of system fll8p . The case of two parabolic and two hyperbolic sectors 
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If polynomial Fi{u) has two real roots which are not equal to zero, it follows from 
(CI), (C2) and Proposition that F2{v) has three real roots and we should consider 
F2{v) instead of Fi{u) (analogous to the case of three real roots of Fi{u)). 

Polynomial Fi[u) can have two real roots, and one of them u = 0. This necessarily 
implies (together with Proposition ([6])) that F2{v) also has two real roots, one of them 
■0 = 0. So we have three roots {u corresponds to v = 1/u which of the same type). 
Simple arguing shows that these three roots can be (1) three saddles (six hyperbolic 
sectors); (2) two nodes and a saddle (two elliptic and two parabolic sectors); (3) a node 
and two saddles (two parabolic and two hyperbolic sectors). 

If Fi{u) has no real root it implies that F2{v) has the only root -0 = and this case 
is similar to the case of one real root of Fi{u): two hyperbolic or two elliptic sectors 
depending on the type of {0,v = 0). 

Other cases are not possible in system (fTHl) with n = 2 that satisfies (CI) and (C2), 
which completes the proof. □ 

If we specify some additional conditions on (|T71) . we can say more about asymptotics 
to the lines y = and a; = of vector field ( ITTl) . Namely, if 

(C3) P„(0, y) = 0^ P*(0, y) = 0, Q„(x, 0) = ^ Q*(x, 0) = 

hold, we can prove 

Proposition 8. Let y = consist of orbits of non-degenerate vector field f|T7j) satisfying 
(C3), and [3 = qn,i/Pn,i > 1- Then ffT7|) has a family of O- orbits y = c|x|'^(l + o(l)), 
where c is an arbitrary constant. 

Proof. If y = consists of orbits of vector field 0171) then m = is a root of Fi{u) = 
—po^n+iu"'~^^ + ■ ■ ■ + {,(ln,i — Pn,i)u + Qn+ifi- It implies that qn+i,o = 0, and y = is a factor of 
Qn{x, y), from which Q*{x, 0) = follows (due to (C3)). From (CI) it follows that pnA 7^ 
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0. According to Proposition [3] the hyperbolic equihbrium point (x = 0, n = 0) is a node 
because we assumed that qn,i/Pn,i > 1 which means F^'(0)P„(1, 0) =Pn,i(?n,i ~ Pn,i) > 0. 
There is a family of characteristic orbits of the form u = c|x|^~^(l + o(l)) where c is an 
arbitrary constant. Returning to variables x, y = ux completes the proof. □ 

Using the second blow-up transformation (120|) we obtain 

Proposition 9. Let x = consist of orbits of non- degenerate vector field ffTTl) . and 
(3 = Pi^n/<li,n > 1- Then vector field ( ITTl) has a family of O- orbits x = c|?/|^(l + o(l)), 
where c is an arbitrary constant. 

The proof of Proposition [9] is similar to the proof of Proposition [8l 

Proof of Theorem The proof follows from Theorem [7] and Propositions M and M □ 

General remarks. Vector field (fTTI) is a particular case of so-called vector field with a 
fixed Newton diagram (e.g., pES]). The conditions of existence of orbits with characteristic 
directions for vector fields with an arbitrary Newton diagram are also known [321 ESI EZ] • 
The following facts are stated for a homogeneous case but can be carried onto more general 
situations almost literally. 
Let 

d d 

ox oy 

be a homogeneous vector field having characteristic directions. The following theorem 
holds. 

Theorem 8. There exists a neighborhood Vt of singular point 0(0,0), such that vector 
field ( IT7|) is topologically equivalent to vector field Jn inQ. 

In case when J„ has no orbits with characteristic directions point 0(0,0) is mon- 
odromic (this can be the case only if n is odd). If we assume that 

(04) Polynomial F{x, y) has no multiple complex factors 

holds, then the cases of focus or center can be identified. 
Denote 

I = > res^— — 

^ F(l,u) 

Im2>0 ^ ' ^ 

the generalized first Lyapunov value. 

Theorem 9. Let 0(0,0) be a monodromic singular point of non- degenerate vector field 
Jn satisfying (04). Then 0(0,0) is a focus if I ^ 0. This focus is attracting if I < and 
repelling if I > 0. 

Thus, in principle, it is possible to analyze non-hyperbolic equilibrium 0(0,0) of sys- 
tem (ITSl) for an arbitrary n. 
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